Precessing supermassive black hole binaries and dark energy measurements with LISA 
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Spin induced precessional modulations of gravitational wave signals from supermassive black hole 
binaries can improve the estimation of luminosity distance to the source by space based gravitational 
wave missions like the Laser Interferometer Space Antenna (LISA). We study how this impacts the 
ablity of LISA to do cosmology, specifically, to measure the dark energy equation of state (EOS) 
parameter w. Using the ACDM model of cosmology, we show that observations of precessing binaries 
with mass ratio 10:1 by LISA, combined with a redshift measurement, can improve the determination 
of w up to an order of magnitude with respect to the non precessing case depending on the total 
mass and the redshift. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 



When the proposed orbiting Laser Interferometer 
Space Antenna (LISA) detects an inspiralling compact 
binary system, it can not only localize the source on the 
sky but can also measure its luminosity distance indepen- 
dent of astronomical distance ladder calibrations. If an 
electromagnetic (EM) counterpart associated with this 
GW event provides the redshift to the source, then the 
combination of these observations can have profound cos- 



mological implications 
of Hubble's constant [l|, 



L| .such asprecise determinations 
., _I EL EL HL H 01 and measurements 
of the dark energy equation of state parameter w [1, Q . 

Spin effects may also help to improve cosmological 
measurements using gravitational waves. If the com- 
pact binary components have non-aligned spins, then the 
modulations induced by precession [l(| [ll| can break de- 
generacies between various parameters being estimated 
and improve accuracy, especially of the sky location and 
luminosity distance [12, 0, [H| ■ Recently we have devel- 
oped a code to carry out parameter estimation for pre- 
cessing inspiralling massive binary black holes, using the 
Fisher matrix formalism, in order to consider the im- 
pact of spin precession on LISA's ability to distinguish a 
general class of massive theories of graviton from general 
relativity [l5j| (see also a similar work by [111). In this re- 
port, we use a variant of this code to study the cosmolog- 
ical implications of the improved distance measurements 
possible with spinning massive black hole binaries. 

Vecchio [12J first pointed out the possible improve- 
ments in precision provided by precessions induced by a 
subset of spin-orbit couplings. Lang and Hughes [H, [3| 
generalized this to the full panoply of spin orbit as well 
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FIG. 1: la errors in the dark energy EOS parameter w as 
a function of the total binary black hole mass at redshifts of 
0.5, 1 and 1.5. Binaries are all assumed to have precession 
and contain black holes of mass (1,10) x 10 4 , (1,10) x 10 5 , 
(5,50) x 10 5 , (1, 10) x 10 6 , (5,50) x 10 6 M Q . The data points 
are medians of 10 4 runs for each mass and redshift, and ignore 
the effect of weak lensing. 



as spin-spin effects. Ref. [l3( focused on the improve- 
ment in the estimation of masses and spins of the binary 
and briefly discussed improvements in distance measure- 
ments. Their follow-up paper [Tij showed that precessing 
binaries would offer much better angular localization by 
LISA and discussed how electromagnetic follow-ups could 
be used effectively to identify the host galaxy and obtain 
the redshift (see also Ref. |17|). 

In this paper, we show explicitly that improved dis- 
tance measurements with precessing binaries combined 
with a redshift to the source could lead to precise mea- 
surements of w; the results are summarized in Fig. [T] 
For example, for a binary system of (1 + 10) x 10 6 M Q 
at redshift z = 1.5, the median la error in measuring w, 
over an ensemble of 10 4 binaries distributed randomly in 
spins, orbital orientations, and sky locations is about 2 
percent. 
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II. MODELS AND ASSUMPTIONS 

Our waveform model is described in detail in [l5| . We 
use second-post-Newtonian (2PN) accurate "restricted" 
waveforms (RWF) for binaries on quasi-circular orbits 
and include the precession effects of spin-orbit and spin- 
spin coupling. The stationary phase approximation is 
used for computing the Fourier transform of the signal. 
Though spin contributions at 2.5PN [l8| and nonspin- 
ning terms up to 3.5PN order [l9| are available for the 
RWF, we assume that they induce only corrections to 
the leading effects studied here. However, incorporation 
of spin dependent higher harmonics [111, [2(| could have 
some influence on the results [2l|. We follow Cutler [22j 
in our model of the LISA satellite and its orbital motion. 
The noise characteristics of LISA that we assume and 
use are the same as in Ref. [23j and used by Lang and 
Hughes [l4|- We have taken into account the effect of 
precession on the antenna pattern functions of LISA. 

We use a Fisher matrix analysis to estimate the er- 
rors in estimating the 15 parameters that characterize 
the system: two masses, two dimcnsionless spin magni- 
tudes (which vary from to 1), the time and phase of 
coalescence, four sets of two angles each specifying the 
location of the binary, the initial angular momentum di- 
rection and the two initial spin directions of the binary's 
members (eight in total), and finally the luminosity dis- 
tance. All the angles used are with respect to the solar 
system bary center. In the specific case that the individ- 
ual spin vectors of the black holes arc aligned, (a rather 
optimistic case astrophysically) the two extra sets of an- 
gles (four parameters) for the individual spins are not 
needed since the binaries do not precess [23]. We also 
assume that LISA provides two independent signal out- 
puts with uncorrelated noises. Finally, we assume that 
the sources are observed for one year prior to coalescence. 

For a given choice of the physical masses of the two 
black holes and of the redshift or luminosity distance, 
we distribute 10 4 sources randomly in the sky with ran- 
dom values of the remaining 10 parameters (we choose 
coalescence time and phase to be one year and zero, re- 
spectively, in all cases). For each of the realizations, we 
solve numerically the precessing equations during the in- 
spiral phase of the system, compute the output signals 
h 1 , h 11 , and the Fisher information matrix defined as, 



1 89- 1 d9 b 



+ 



dh 11 



dh u \ 

de b J 



where the inner product is, 



(ftlWE4Re l dfi s ni f) 



(2.1) 



(2.2) 



where h(f) denotes the Fourier transform of the gravita- 
tional waveform h(t, 9 a ), star denotes complex conjugate, 
S n (f) is the noise spectral density of the detector, and 
dh/d9 a denotes the partial derivative with respect to the 
parameter a being estimated. The superscripts / and 




FIG. 2: Distribution of the lcr errors in luminosity distance 
measurement based on 10 4 random realizations. The mass of 
the binary system is (10 5 + 10 6 )Mq at a redshift of z = 1. 
Solid (red) histogram is for precessing systems; dashed (blue) 
histogram is for systems with spins aligned. 



II denote the two LISA outputs. The signal-to-noise 
ratio (SNR) for a given signal h(t,9 a ) is then given by, 
p[h] = (h^) 1 / 2 . We then invert the Fisher matrix to 
obtain the covariance matrix E ah , and the correspond- 
ing root mean square errors from the square roots of the 
diagonal entries, as follows, 



A9 a = 



s = r~ 



(2.3) 



We focus on systems with a mass ratio of 10 : 1 since 
these are astrophysically interesting and exhibit stronger 
effects of spin precession. 

Fig. [5] shows the histogram of the relative errors 
in the luminosity distance for a system with masses 
(10 5 ,10 6 )M Q at a redshift of z = 1. Precession of the 
binary has a significant impact on the distribution as the 
peak of the distribution shifts to the left (improves) al- 
most by an order of magnitude with respect to the case 
where all systems in the population are nonprecessing. 
This is in reasonable agreement with the results of Lang 
and Hughes [l4| (see, e.g., figure 7 of their paper). 

However, to do cosmology, we need the redshift of an 
electromagnetic signal associated with an observed LISA 
event or with its host galaxy. To date there is no clear un- 
derstanding of the astrophysical mechanisms that could 
produce an electromagnetic afterglow (or a precursor) to 
a binary black hole merger, although a number of pos- 
sibilities have been discussed (24|. Whether LISA can 
localize the source on the sky so that extensive electro- 
magnetic follow up missions can be launched has been 
addressed in many recent works. Of particular interest 
here are the estimates of Lang and Hughes [l4[ , showing 
that typically for a z = 1 source of total mass 10 6 M Q , the 
angular resolution taking precession into account is about 
0.3 — 3 square degrees one day prior to merger (see Table 
5 of [3] for example). Kocsis et al [25| argued that find- 
ing EM counterparts associated with LISA sources will 
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be difficult, but may be achievable with the advent of 
various wide field telescopes which would be operational 
by the time LISA flies (see Table I of Refj2f| for details). 
Further, some authors have argued [1, l24l] that the num- 
ber of candidate galaxies associated with a LISA event, 
can be reduced by 2-3 orders of magnitude by looking for 
the source within a 3D error volume, combining the an- 
gular resolution of LISA and the approximate luminosity 
distance LISA would provide in advance. For this reason 
(rather optimistically) we have not rejected any distance 
estimate in the 10 4 realizations based on the size of the 
angular error box or the detectability of an EM counter- 
part. 

To translate distance errors into errors in the dark 
energy EOS parameter w, we work with the standard 
cosmological model with a flat universe and the nomi- 
nal parameters Ho = 75 kms" 1 Mpc -1 , matter density 
Qm = 0.25, dark energy density Q\ = 0.75 and w = — 1. 
The luminosity distance of a source at a redshift z is 
given by, 



D, 



(1 + *) 



dz' 
o W)' 



(2.4) 



where H(z) is given by 



H(z) = H oy /n M (l + zf + fi A (l + z)3( 1 +-) . (2.5) 

In order to illustrate directly the effect of precession, 
we assume that errors in Hq, Qm, and z are negligible, 
so that the error in w is related directly to the error in 
D L by HH 



Aw = D L 



dD, 



dw 



-1 AD, 



(2.6) 



and 



The value of dD^/dw can be calculated using (|2.4[) 
(|2.5p for the different values of the redshift used. 

Finally, we have neglected the effect of weak lensing 
on the estimation of the luminosity distance. This is 
a serious caveat, especially at redshifts above 1. How- 
ever, mechanisms have been proposed which might help 
to reduce the weak lensing error in the future. These in- 
clude "corrected lenses" [231, the accurate measurement 
of the weak lensing power spectrum with radio observa- 
tions p8j . and combining optical and infra-red observa- 
tions of foreground galaxies [29| . See [3(| for a recent 
discussion. 



III. RESULTS AND DISCUSSION 

The inclusion of spin precession and its attendant mod- 
ulation of the GW waveform has a significant impact on 



the accuracy of measurements of the dark energy param- 
eter, which in turn is due to improved estimation of the 
luminosity distance, as shown in Fig. [21 Figure [1] depicts 
how the estimation of w varies with the total mass of the 
binary. We show the median errors from various runs 
corresponding to different masses and redshifts. It is in- 
teresting to note that for masses between 10 5 - 1O 7 M , 
Aw < 3% for redshifts up to 1.5. The left panel of Fig. [3] 
shows that the errors in the measurement of w go down 
by a factor of approximately 10 when the spins of the 
binaries are not aligned compared to the aligned, non- 
precessing case. Note that all our histograms, except 
the cumulative plots in the insets, are unnormalized, so 
that the histograms show the number of realizations in 
each bin. The dependence of the errors in the redshift 
z is shown in the right panel of Fig. [3] As is expected 
the errors get worse with increasing redshift (distance), 
but the degradation of errors is not dramatic. Even for 
z = 1.5, a significant fraction of the 10 4 randomly dis- 
tributed sources permit measurement of w to better than 
5% accuracy. However, for larger redshifts, more galax- 
ies will be found in the LISA error volume and it will be 
more difficult to obtain a redshift. 

Van den Broeck ct al. [3l[ have recently revisited the 
measurement accuracy of the dark energy EOS parame- 
ter w, using a GW signal without precession but includ- 
ing higher harmonics. Though higher harmonics and pre- 
cession are completely different effects, it is interesting to 
note that both improve the estimation of w significantly 
and to comparable orders. This further motivates the at- 
tempt to incorporate the effect of higher harmonics in the 
presence of spin precession (see [2l[ for a recent analysis). 
We have ignored any possible effect of orbital eccentricity 
in the waveforms (see, eg. [32[ and references therein). 
These are topics for future consideration. 
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FIG. 3: Distributions of la errors in w for various values of redshifts and spin configurations. Left panel: Distribution of 
the errors in the dark energy EOS parameter w for a binary system of masses 10 5 + 10 6 M@ at redshift z = 1, and for the 
precessing and non-precessing cases. The inset frame shows a zoom of the corresponding cumulative histograms. Right panel: 
Distribution of errors in w for the same binary system with masses 10 s + 10 6 Mq in the precessing case for three different 
redshifts. The inset frame shows a zoom of the corresponding cumulative histograms. 
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